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Abstract 

Using 3.07 pb~^ of data collected in the energy range 0.60-0.97 
GeV by CMD-2, about 150 events of the process e'^e" — > vr+Tr^Tr+vr" 
have been selected. The energy dependence of the cross section agrees 
with the assumption of the ai(1260)7r intermediate state which is dom- 
inant above 1 GeV. For the first time tt~^ Tr~ ir'^ Tr~ events are observed 
at the p meson energy. Under the assumption that all these events 
come from the p meson decay, the value of the cross section at the p 
meson peak corresponds to the following decay width: 

r(/>° vr+vr^vr+vr") = (2.8 ± 1.4 ± 0.5) keV 

or to the branching ratio 

B{p^ vr+vr^vr+TT") = (1.8 ± 0.9 ± 0.3) • 10~^ 

1 Introduction 

Production of four pions in e~^e~ annihilation is now well studied in the 
cm. energy range 1.05 to 2.5 GeV (see [ |l| and references therein). Much 
less is known about this process at energies below the meson mass. The 
first attempts to observe it resulted in detection of single events above 0.96 
GeV by Orsay groups [ |], || while scans of the energy range 0.64 to 1.00 GeV 
by Novosibirsk groups placed upper limits only [ |, |^. These observations 
confirmed the predictions that the cross section is very small at these energies 

The situation was improved in 1992 when the upgraded high luminosity 
collider VEPP-2M resumed its operation at the Budker Institute of Nuclear 
Physics in Novosibirsk [ With two modern detectors CMD-2 [ ^ and SND 
[ |TD[ various high precision measurements in the cm. energy range from the 
threshold of hadron production to 1.4 GeV became possible. CMD-2 has 
recently published results of their analysis of the four pion production above 
the (j) meson The large data sample allowed a precise measurement of the 
energy dependence of the cross section and revealed the intermediate states 
through which final pions are produced. It was shown that while the reaction 
e^e~ 71^71^271^ proceeds through the utt^ and 01(1260)7? intermediate 
states, it is the latter which dominates in the vr+Tr^Tr+Tr^ channel. 

In this paper, we extend the analysis of the process e^e~ —>■ ti^ti^ti^tt^ to 
the cm. energy range from 0.60 to 0.97 GeV. The high integrated luminosity 
allowed observation of 153 ti^ti^ti^tt^ events and the first measurement of 
the corresponding cross section at a level as low as several tenths of pb. 
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2 Experiment and data analysis 



The experiment was performed in spring 1998 when the energy range from 
0.97 down to 0.60 GeV was scanned with a typical step of 10 MeV outside 
the p and uj resonances and 1 MeV near the peak of the resonances. To 
decrease statistical errors, in the final analysis the information from several 
low energy points was combined. The two first columns of Table 1 present 
the corresponding energy points (intervals) and integrated luminosities. The 
total integrated luminosity was 3.07 pb^^. 

The general purpose detector CMD-2 has been described in detail else- 
where [^. It consists of a drift chamber and proportional Z-chamber, both 
used for the trigger, and both inside a thin (0.4 Xq) superconducting solenoid 
with a field of 1 T. 

The barrel calorimeter placed outside the solenoid consists of 892 Csl 
crystals of 6 x 6 x 15 cm^ size and covers polar angles from 46° to 134°. The 
energy resolution for photons is about 9% in the energy range from 50 to 600 
MeV. 

The end-cap calorimeter placed inside the solenoid consists of 680 BGO 
crystals of 2.5 x 2.5 x 15 cm^ size and covers forward-backward polar angles 
from 16° to 49° and from 131° to 164°. The energy and angular resolution 
are equal to (Je/E = A.6%/^E{GeV) and cx^^e = 2 ■ IQ-^ / ^fE{GeV) radians 
respectively. 

The luminosity was determined from the detected e+e~ — > e+e~ events [ 

For the analysis of the reaction e+e^ —>■ 7r+7r^7r+7r^ events with four 
charged tracks were selected. The tracks were required to originate from the 
interaction region: 

• the impact parameter of the tracks r^m is less than 0.3 cm 

• the vertex coordinate along the beam within ±10 cm. 

To have good reconstruction efficiency, tracks were also required to cross at 
least two layers of the drift chamber: | cos6'| < 0.8. 

For the selected events the kinematic fit was performed assuming that all 
tracks are pions and under the constraint that the sum of the 3-momenta 
YjiVi = 0- Further analysis was performed using the normalized "apparent 
energy" : 

Ei + ml 



3 




0.8 



0.5 



1.5 



2.5 3 

T mm 



Figure 1: Distribution in the normalized apparent energy versus the mini- 
mum space angle between tracks with the opposite charges 



Figure |I] shows the distribution of Eapp versus the minimum space angle 
between the tracks with opposite charges ipmin- A narrow band with Eapp ~ 1 
corresponding to 'k^tt~t[^tt~ events is clearly observed in the region "a" . The 
region "b" is populated by events from the processes: 

e^e~ — > 7r'''7r~7r°, vr" — > e^e~'y (1) 

and 

e"^e^ e"'"e~7, (2) 

e^e" 7r^7r~7 (3) 

with the subsequent conversion of the photon into a e^e^-pair at the beam 
pipe. 

Events of the process e'^e^ tt+tt^, where the products of the nuclear 
interaction of the pions scatter back into the drift chamber and induce two 
"extra" tracks, fall in the region "c" . To suppress the background from this 
process, the condition Pnorm < 0.8 was imposed where the parameter of the 
normalized maximum track momentum was defined as 

maxi \pi\ 



P 

^ norm 
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To reject the background from the reaction (|I]) in the p/uj meson region 
where the number of observed events is small, the e/vr-separation procedure 
from [ Q was employed. The idea of the method is the following. Using 
reconstructed events of the processes e^e~ —>■ 7r"'"7r~7r° and e~^e~ —>■ e~^e~'~f in 
the (j) meson energy range, the distributions in the parameter Ecsi/\p\ were 
studied for both particle types (e/vr) and signs (+/—). Here Ecsi is the en- 
ergy deposition in the Csl calorimeter and p is the momentum of the particle 
with a track matching the cluster in Csl. The probability density functions 
fn+ ) f-iT- , fe+ and fe- for the parameter Ecsi/IpI were determined for 50 MeV 
bins of particle momenta (0-50, 50-100 etc.). Then the probabilities were 
defined: 



fn± + fe± /vri + /e± 

for a track to be respectively pion or electron. In the same way the probability 
was defined: 

w 



for the pair of tracks to be pions. Figure ^ shows the distribution in the pa- 
rameter Wt^+tj- for the events from the processes e+e~ e+e~7 (histogram) 
and e+e~ tt+tt^tt^ (points with errors) in the (f) meson region. As is clear 
from this Figure, the parameter efficiently separates events with e+e~ 

pairs from those with tt+tt" pairs. In our analysis we calculated the param- 
eter W-n+T,- for two tracks with the smallest angle between them. Figure ||b 
demonstrates the distribution in this parameter for the events from the re- 
gion "a" (in Figure 1^) which is dominated by 7r+7r^7r+7r~ events. We required 
Wt,+t,- > 0.5 to suppress the background from the process (|I]) in the energy 
range 0.6 to 0.8 GeV. The efficiency of this selection criterion was determined 
using the 7r"'"7r~7r"'"7r~ events in the energy range 0.81 to 0.97 GeV. 

Finally, the selection of events with Sapp < 1.05 and ilJmin > 0.3 provides a 
practically pure sample of 153 7i^n^7i^7i~ events. The remaining background 
from the reactions (P, (^ and (|^) was estimated using the complete Monte 



Carlo simulation (MC) of the CMD-2 detector [ |T3[ and appeared to be 
negligible. For example, in the energy range 0.6 to 0.8 GeV the expected 
number of events from the process (|l|) is less than 0.2 to be compared to four 
events observed. 

Figure ^ shows the distribution of the selected events in Sapp- In this 
Figure black squares with error bars represent the experimental data while 
the shaded histogram is simulation of the process e~^e~ 7r"'"7r~7r"'"7r~. Good 
resolution in the parameter Sapp allows efficient separation even of the exotic 
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Figure 2: e/vr separation probability, a — distribution in for the 

events from the processes e^e~ —* e^e~7 (histogram) and e'^e^ tt+tt^tt'^ 
(points with errors) in the meson energy range; b — distribution in Wt^+t^- 
calculated for two tracks with the smallest angle between them for the 
7r+7r^7r+7r~ events (region "a"). The arrow shows the restriction imposed 
upon the parameter M/V+tt- 



final state 7r"'"7r~/x"'"/i~. The distribution for it is shown by the transparent 
histogram obtained by the simulation of the process e'^e" 7r+7r~/i"'"/i~ 
assuming the phase space kinematics of the final particles. 



3 Determination of cross sections 

The cross section of the process e+e~ 7r+7r^7r+7r^ was calculated at 
each energy using the formula: 

Li Ei (1 + bi) 

where iVj is the number of selected 7r+7r^7r+7r^ events, Li is the integrated 
luminosity, Ei is the detection efficiency, and 5j is the radiative correction at 
the z-th energy point. 

The detection efficiency was determined from MC assuming four differ- 
ent quasi-twobody production mechanisms ai(1260)7r, pa, 7r(1300)7r and 
02(1320)71. In contrast to our observations at the energy above 1.05 GeV [ 
H, in the energy range studied in this experiment the effects of interference 
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Figure 3: Distribution in Sapp- Separation of the processes e+e 
7r"'"'7r~7r+7r^ and e"'"e~ 7T~^'n'~ fi~^ fi~ 



between amplitudes differing by permutations of identical pions make the dif- 
ference between the models smaller. Comparison of various invariant mass 
and angular distributions predicted by the model with experiment showed 
that 7r(1300)7r and a2(1320)7r intermediate states can not well describe the 
data whereas ai(1260)7r and pa are almost indistinguishable in this energy 
range. For the final calculations the ai(1260)7r production mechanism which 
clearly dominates at higher energy [ |I| has been chosen. The detection effi- 
ciency grows with energy varying from 17 to 28 % in the energy range 0.75 
to 1.0 GeV. It falls quickly below 0.7 GeV since pion momenta become too 
small to reach the Z-chamber and produce a trigger. 

Radiative corrections were calculated according to [ [T4|. 



Table Q gives the summary of the cross section calculations. Since the 
number of selected events below 0.8 GeV is small, the data from separate 
points were combined. The values of the cross sections or upper limits were 
obtained from the following formula: 

^ " SL,e,(l + 5,) 

which takes into account the variation of the luminosity, detection efficiency 
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Table 1: Summary of the cross section calculations 



E 

GeV 


nb-i 




£ 


5 


0", 

nb 


0.6-0.66 


122.3 








029 


-0 


147 


< 0.76 at 90% CL 


0.69-0.72 


103.2 








115 


-0 


157 


< 0.23 at 90% CL 


0.75-0.78 


372.8 


2 





171 


-0 


134 


0.036 ± 0.026 


0.781-0.8 


818.3 


2 





181 


-0 


133 


0.015 ± 0.011 


0.81 


61.3 


2 





243 


-0 


117 


0.15 ± 0.11 


0.82 


112.1 


3 





217 


-0 


114 


0.14 ± 0.08 


0.84 


133.7 


4 





238 


-0 


113 


0.14 ± 0.07 


0.88 


172.1 


9 





258 


-0 


116 


0.23 ± 0.08 


0.92 


292.7 


26 





275 


-0 


121 


0.37 ± 0.07 


0.94 


140.2 


10 





275 


-0 


124 


0.30 ± 0.09 


0.95 


232.0 


18 





284 


-0 


126 


0.31 ± 0.07 


0.958 


256.6 


34 





284 


-0 


127 


0.53 ± 0.09 


0.97 


256.4 


43 





276 


-0 


130 


0.70 ± 0.11 



and radiative correction over the energy interval and gives basically the cor- 
rect average cross section. 

Figure ^ shows the energy dependence of the cross section below 1.05 
GeV. Only statistical errors are shown. For illustration, also demonstrated 
in the Figure are results from other measurements [ ^, ^, ^, ^ |15[. The values 
of the cross section obtained in this work are consistent with them and match 
the measurements of CMD-2 above the meson. 

The systematic uncertainty comes from the following sources: 

• selection criteria - 10 % 

• event reconstruction - 5 % 

• detection efficiency dependence on production mechanism - 3% 

• luminosity determination - 1.5 % 

• radiative corrections - 1 %. 

The overall systematic uncertainty was estimated to be ~ 12%. 
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Figure 4: Cross section of the process e~^e h^ti h^ti below 1.05 GeV 
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Figure 5: Cross section of the process e^e" 7r+7r~7r+7r^ in the cm. energy 
range 0.65-2 GeV. Comparison of our data to the data of other groups. The 
dashed and dotted curves show separately the contribution of the p and p' 
mesons respectively while the solid one is their sum taking into account the 
interference between them 
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4 Discussion 



The shaded area in Figure |^ corresponds to the extrapolation of the en- 
ergy dependence of the cross section from the energy region above 1.05 GeV 
[ m . The calculation assumed that the cross section behaviour is determined 
by two interfering resonances - p(770) and its excitation (p') decaying into 
the final four pion state via the ai(1260)7r intermediate mechanism. The 
coupling constant Qpanr was fixed from the width of the ai(1260) decay to 
p-K. Parameters of the p' were determined from the fit of the high energy 
data of CMD-2 [ |^ and DM2 [ |l^]. The central curve corresponds to the 
ai(1260) width of 600 MeV optimal in our analysis whereas the upper and 
lower curves are obtained for the widths of 800 and 400 MeV respectively. 

It can be seen that the energy dependence of the present data is consistent 
with the assumption of the ai(1260)7r dominance earlier established at higher 
energies [ |l| . 

Figure ^ shows the energy dependence of the cross section o"(e+e~ 
7r+7r~7r+7r~) in the energy range 0.65-2 GeV. In addition to our results also 
shown are the cross sections determined by various groups in Frascati, Orsay 
and Novosibirsk [ |, |, ||, |l|, |17l |18|, |T|, |20|, [|. The dashed and dotted 



curves show separately the contribution of the p and p' mesons respectively 
while the solid one is their sum taking into account the interference between 
them. 

It can be seen that at low energies the p' contribution is much less than 
that of the p (about 2% at 0.8 GeV). However, it quickly grows with energy, 
reaching the level of 15% already at 1 GeV. 

We do not observe the energy dependence characteristic of the resonance 
with the Breit-Wigner fall of the cross section at the right slope of the p 
meson. Instead, a shoulder can be seen since the resonance behaviour is 
compensated with the fast growth of the phase space of four pions. To 
compare the obtained experimental data to existing theoretical predictions, 
it is useful to evaluate the cross section of the process at the cm. energy 
corresponding to the p meson mass. To this end the energy range 0.75-0.8 
GeV was considered. Using the information about the integrated luminosities 
in this range and following the same procedure of averaging, the cross section 
was determined as above. The value of the cross section at the p meson mass 
is (0.020 ± 0.010 ± 0.003) nb. This value can be conventionally presented in 
terms of the width of the p meson decay to 7r+7r^7r+7r^. 

Using the expression for the cross section of the resonance production at 
its peak 

^0 = ^2 (6) 
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and the values of M and B^e from [ |2T]] one obtains the branching ratio 

5(p° ^ 7r+7r"7r+7r") = (1.8 ± 0.9 ± 0.3) ■ 10"^ 

or for the width 

r(p° ^ 7r+7r"7r+7r~) = (2.8 ± 1.4 ± 0.5) keV. (7) 

to be compared with the best previous measurement in which an upper limit 
< 30 keV was placed [ g. 

Attempts to describe the four pion production are known since rather 
long time. Earlier theoretical models devoted to the calculation of the cross 
section of e^e^ 11^71^71^71^ , were based upon the natural assumption of 
the quasi-twobody enhancement of the p tail [ 0]. In most of them the 
predictions are given for the energy range above 1 GeV, while the cross section 
falls very rapidly when the energy decreases below this energy. For example, 
in [ 1^ the cross section of 7i~^7i~7i~^7i~ production is calculated assuming 
ai(1260)7r, a2(1320)7r and pa intermediate states. The ai(1260)7r is predicted 
to dominate above 1.05 GeV, while only the pa mechanism survives at lower 
energies. The absolute magnitude of the cross section due to the above 
mechanisms is about 1 pb at 1.05 GeV and much smaller at lower energies, 
in obvious conflict with our observations. This model also fails to reproduce 
the cross section at higher energies typically predicting an order of magnitude 
smaller values than observed. In most of these papers interference effects are 
not taken into account leading to wrong relations between the probabilities of 
the final states differing in the number of charged and neutral pions. Another 
serious defect of all calculations was the lack of the p' which was unknown 
at that time and which plays a rather important role at the energies above 
1 GeV. Finally, the expressions for the matrix elements contain form factors 
with unknown energy dependence and their wrong choice can seriously affect 
the magnitude of the cross section. 

While theory fails to explain phenomena far from the threshold, one could 
expect that in the region of energies close to the p meson peak there would be 
more successful calculations. Various theoretical predictions [ |22|, |23|, 
26| for the value of r(p° 7r+7r^7r+7r^) are compared to our result in Table |^. 



Most of them use the low energy effective Lagrangian approach. An estimate 



based on the broken SUexOs quark model [ 22] did not incorporate chiral 



symmetry and was ruled out by the measurement in Ref . [ m . Later attempts 



to include the p meson in the Lagrangian [ ^ |2^ also failed to completely 
respect the chiral symmetry as noted by the authors of Ref. [ |2^ in which 
various effective Lagrangians for vr and p (and in some cases ai(1260)) were 
suggested. The predictions of Ref. [ ^ are an order of magnitude smaller 
than these in \ ESI 12^ and do not contradict our measurement within errors. 
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Table 2: The decay width r(p° tt+tt tt+tt ) 



Work 


Tp''— >7r+7r 7r+7r ) keV 


[ 22| 


172 


[ 23| 


7.5 ± 0.8 

25 ± 3 
60 ± 7 


[ 


24| 


16 ± 1 


[ 


25] 


0.59-1.03 


[ 


26| 


0.89 


This experiment 


2.8 ±1.4 ±0.5 



The authors of a recent pubhcation [ |2^ also describe the p ^ Air decay 
using the chiral Lagrangian (not including the ai(1260) meson) and come 
to the value of T{p^ —>■ 7r+7r^7r+7r^) close to that in [ In contrast to [ 
p5| , they consider the cross section rapidly varying with energy as a more 
adequate characteristics of the process than the decay width and calculate the 
value of the cross section up to 0.9 GeV. Above 0.8 GeV their predictions for 
the cross section are well below the measured values. Apparently, the chiral 
Lagrangian approach fails to describe the dynamics of four-pion production 
already at these modest energies. The authors of Ref. [ ^ ascribe the 
observed discrepancy to the fact that their model neglects higher derivatives 
and chiral loops. However, from the phenomenological point of view it is clear 
that the existence of the p' meson is essential for understanding the observed 
pattern while chiral theories fail to include excited states. In general, chiral 
theory seems to produce valid predictions at very low energies only. 

As already noted above, thorough analysis of the four pion production 
above 1.05 GeV allowed to establish the ai(1260)7r dominance in e~^e~ anni- 
hilation [ |l| and successfully describe a four pion decay mode of the r lepton 
[ E^, |28|. We could also observe that the extrapolation of the model pre- 



dictions to lower energies is in good agreement with the experiment even at 
energies as low as 0.8 GeV. 

Better understanding of this energy range could come after a bigger data 
sample is collected as well as determination of the cross section of the reaction 
e~^e~ — > 7r~^iT'7r^7T^ and its combined analysis together with the reaction 
studied in this work. 
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5 Conclusions 



For the first time the cross section of the process e^e^ — ^ 7r+7r^7r+7r^ has 
been measured in the cm. energy range 0.60-0.97 GeV in which 153 events 
are observed. The cross section falls rapidly with the energy decrease and 
does not contradict the model in which final pions are produced via p(770) 
and p' mesons decaying into ai(1260)7r. 

The cross section value at the peak of the p meson corresponds to the 
following value of its decay width 

r(p° ^ Tr+TT-TT+TT-) = (2.8 ± 1.4 ± 0.5) kcV 

and its branching ratio 

B(p ^ TT + TT^TT + Tr") = (1.8 ± 0.9 ± 0.3) ■ 10"^ 

The measurement of the R value at the energies below 1 GeV known to give 



the dominant contribution to the hadronic part of the (g-2)^ [ ^ is one of 
the main goals of CMD-2. Results of this work show that we can measure 
small cross sections even at the level of one tenth of nb and can thus provide 
the high precision determination of R. 
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